Conformally coupled dark matter. 



Mark Israelit Q, 



Department of Physics, University of Konstanz, PF 5560 M678, D-78434 Konstanz, Germany. On 
from: Department of Physics, University of Haifa, Oranim, Tivon 36006, Israel. 

Abstract 

Dark matter is obtained from a scalar field coupled conformally to gravitation; the scalar being a 
relict of Dirac's gauge function. This conformally coupled dark matter includes a gas of very light 
(m » 2.25 x 1(T 34 eV) neutral bosons having spin 0, as well as a time-dependent global scalar field, 
both pervading all of the cosmic space. The time-development of this dark matter in the expanding 
F-R-W universe is investigated, and an acceptable cosmological behaviour is obtained. 

Keywords: General Relativity; Cosmology; Dark Matter. 



1 c-mail: Marc.Israelit@uni-konstanz.de, and: israelit@physics.technion.ac.il, permanent e-address 



1 Introduction. 



There is a great deal of dark matter in the universe, the dark matter having a much greater mean 
density then visible matter, (Trimble 1987), (Turner 1991), (Tremain 1992). It is believed (Tyson 1992), 
(Priester, Hoell, & Blome 1995) that the density of dark matter exceeds that of luminous by a factor 
of about 10. At present there is a number of popular non-baryonic dark matter candidates: axions, 
neutrino, gravitino etc. etc. However, these candidates suffer from a number of shortcomings: their 
evolutionary behaviour in the expanding universe is not clearly known, their origin is not related to the 
geometry of the universe, they are inserted into the space-time framework from the outside. 

Some years ago it was suggested (Israelit, & Rosen 1992) that at least a part of the dark matter 
in the universe consists of massive bosons of spin 1. These particles, named weylons were created by 
a local Weylian vector field. The idea of such a field was based on the Weyl geometry (Weyl 1919) as 
modified by Dirac (1973). It was also shown (Israelit, & Rosen 1994) that dark matter consisting of heavy 
weylons (with mass greater then 10 Mev), may be fitted into the standard cosmological model (Weinberg 
1972), as well as into the singularity-free prematter cosmological model (Israelit, & Rosen 1989). Later 
a global dark matter effect was obtained (Israelit, & Rosen 1995) from a time-dependent scalar field, 
the scalar being identical with the gauge function introduced by Dirac (1973) in his modification of the 
Weyl geometry. The above mentioned dark matter forms are derived from the Weyl geometry. In this 
paper an additional geometrically based generator of cosmic dark matter is presented. 

The purpose of the present work is to consider the possibility that a scalar field, coupled conformally 
to gravity, can generate cosmic dark matter. The behavior of this dark matter in the expanding universe 
will be investigated. It will also be shown that unifying the Einstein-Hilbert action with the (reduced) 
action of the Weyl-Dirac theory, one obtains the action of a scalar field coupled conformally to gravitation. 
In this procedure the coupled scalar field appears as a relict of the Dirac gauge function. 

In the present work general relativistic (geometrical) units are used, so that G = 1, c = 1, the 
Hubble constant is given in cm -1 , and the energy density of matter in cm~ 2 . However, in order to 
discuss physical results we turn sometimes to conventional units, noting this fact in the text. The 
conversion rules may be found in (Synge 1966). 



2 The Weyl-Dirac theory and conformally coupling. 

A detailed description of the Weyl-Dirac theory and a discussion of its physical aspects may be found 
in Rosen's work (1982). A brief discussion of the theory in view of the dark matter problem is given in 
(Israelit, & Rosen 1992). In this section we consider the relation between the Weyl-Dirac theory and 
a scalar field coupled conformally to gravity. The notations of (Rosen 1982), and of (Israelit, & Rosen 
1992) will be used. 

The field equations of the Weyl-Dirac theory may be derived from a variational principle (Dirac 1973), 
(Rosen 1982) 

SId = 0, (1) 

with the action 

Id = J M-ff) 1/2 <*V (2) 
and the Lagrangian density Ld given by 

L D = W^W^-0 2 RZ + k0p0* + (k-6)(20w*0 ll +0 2 w ll w») + 2A0 i . (3) 

In expression (^), w tJ/ (x x ) is the connection vector of the Weylian geometry, the Weylian length curvature 
tensor W^ v is defined as W^ v — w^. v — w v (a comma denoting partial differentiation), 0{x x ) is the 
Dirac gauge function, and /3 M = jlJr . Further, A is the cosmological constant, R G a is the Riemannian 
curvature scalar, and k is an arbitrary parameter. In order to get a geometrically based description 
of gravitation and electrodynamics, Dirac took k = 6. In that case may be treated as the vector 
potential of the electromagnetic field, and W M „ becomes the field tensor. By Dirac's choice one has from 
equation (||) 

L D =W liV W fiU - 2 R <T a +&0 a a + 2K0 i . (4) 
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For cosmological cosiderations we can neglect the Maxwell term in , so that 

L D = ~/3 2 R a a +6(3 a (3 rT + 2A(3' i . (5) 

Inserting (g) into (Q) we obtain the action for a conformal space. Further, combining it with the Einstein- 
Hilbert action for gravitation and matter 

Ig = J(K+L m )(-g)^ 2 d 4 x , (6) 

(L m is the Lagrangian density of matter), we get 

/ = JlK + L m + 8Tr(P°f3 a - l -p 2 R° a + -Ap^i-g^^x . (7) 

This is the general covariant action of a scalar field (3(x x ) coupled conformally to gravitation. In the 
following procedure L m will be considered as that of ordinary (luminous) matter, and it will be assumed 
that L m does not depend on (3. On the other hand, dark matter will be obtained from the conformally 
coupled scalar field (i{x). We also shall consider models without the cosmological constant so that in the 
action (Q) we shall set A = 0. 

3 The field equations. 

Varying in(^) the metric tensor g^ v we obtain the Einstein equation 

R„ u - 1/2 g^R° = -8n{T^ + 9 M „) , (8) 

with the Ricci tensor R^ v , and with the energy- momentum density tensor of ordinary (luminous) 
matter 

87r V = g T ^ = Sgl „, ■ ( 9 ) 

The dark matter is represented in (||) by the energy-momentum density tensor of the /3-field 

60 M „ = 40 ll p v -g lu ,p <r p r -2pp ll . v + 2g^f3(3% - f3 2 (R^ - 1/2 g„ u R°) . (10) 
Further, variation with respect to the field (3 yields the following field equation: 

p% + (l/6)/3< = 0. (11) 
Making use of ([ll]), one obtains from ( |l0| ) 

e- = o , (12) 

so that (||) gives 

R% = 8irTZ . (13) 
Taking into account (|T^) , one can rewrite the field equation ([ll]) as 

0% = -^-(3TZ , (14) 

so that the behavior of the conformally coupled field depends on the amount and state of ordinary matter. 
Finally, making use of the contracted Bianchi identity and of fll]), one obtains from (Il0|) 



e^„ = o , (15) 

and from (||) one has 

?;% = o , (i6) 

so that one obtains separated energy-momentum relations, the first, ( |l5| ) for dark matter, and the second, 
(O) for ordinary matter. 



3 



Let us consider a homogeneous and isotropic spatially closed universe described by the F-R-W line 
element 

ds 2 = dt 2 - R 2 ( — + r 2 dd 2 + r 2 sin 2 d dip 2 J , (17) 

with R(t) being the radius of the universe. We assume that the universe is filled with ordinary cosmic 
matter, having an energy density p(t), and a pressure P(t), as well as with dark matter, created by 
the scalar function f3. For the metric given by ( |l7| ) we can write down the Einstein equation (@) with 
H = 0; v = 0; 

^ + ^ = Y (P + eS) , (is) 

and from eq. ([l3]) we obtain 

R R 2 1 47T, 

R + R 2 + R^-Y {P - 3P) ■ (W) 
In addition we have from the energy relation ( |T^ ) 

P+— (p + P) = 0. (20) 

The equations obtained in this section enable us to investigate the scalar field [3 coupled conformally 
to a F-R-W universe filled with ordinary matter. The generation of dark matter by this /3- field may be 
considered either by means of a local approach, or by a global one. 



4 A local approach. Scalar bosons. 

From the global standpoint, in a homogeneous and isotropic universe (3 depends only on the cosmic time 
t. However, one can assume that there exist small regions of inhomogeneity, in which (3 may depend also 
on spatial coordinates. Introducing local Lorentzian coordinates (cr, x, y, z), one can rewrite ( pr| ) as 

□/3+y(p-3P)/3 = 0. (21) 

Further in small space-time regions one can neglect the dependence of p and P on the cosmic time t, so 
that 

4tt 

— (p - 3P) = k — const. , (22) 
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and ( |21| ) may be rewritten in the form of a Klein-Gordon equation 

D/3 + k 2 (3 = 0. (23) 

From the point of view of quantum mechanics, eq. (|23| ) describes spinless particles, scalar bosons, with 
a mass m and Compton wavelength A c . In conventional units one has 

mc 27T . , 

K =J^ = -> (24) 



and by (123), and (p4) one obtains for the mass in conventional units 



h 

m = — 



47T 

— G{p-iP) 



1/2 

(25) 



with G being Newton's gravitational constant. It will be recalled that p, and P are the energy density 
and pressure of ordinary (luminous) cosmic matter. 

In order to estimate the boson mass at present we will make use of the Hubble constant. Following 
(Peebles 1993) let us write for its present value H N = Hoh, where Hq = 100 km s _1 Mpc^ 1 
(9.78 x 10 9 y)" 1 = 1.08 x 10" 28 cm" 1 (cf. (Weinberg 1972)), and with 0.5 < h< 0.8 (cf.(Priester, 
Hoell,& Blome 1995)) . Adopting the value h ~ 0.5 we get 

H N = 0.5 x 10~ 28 cm" 1 , (26) 
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so that for the critic density we obtain 



Pc 



Hi = 0.3 x l(r 57 cTO- 2 . 



(27) 



To discuss a closed universe one can take for the total ( luminous and dark ) matter density at 
present Ptotai n = 0.35 x 10~ 57 cm- 2 > p c . Further one can assume that the luminous matter 
makes up 0.1 of the total one, so that the following values must be substituted into (Eq): 



P n = , p N = 0.35 x 10~ 58 cm~ 2 
Expressing p N in conventional units one obtains the present value 

m„ re 4 x 10~ 67 = 2.25 x 10~ 34 eV 



(28) 



(29) 



According to(25) in the radiation period the boson mass vanishes, whereas during the dust-dominated 
period it is given by 



m = m. 



R 



1/2 



(30) 



with R N being the present value of the cosmic scale factor. 

Let us consider the energy-momentum density of the dark matter. For small regions, in local 
Lorentzian coordinates, one obtains from ( |Io| ) the following non-vanishing components: 



&0 = ^- c P T P X + j- c Wr X 



(31) 
(32) 



1 



' 2 R R 2 



6 13 1 R + R 2 + R 2 



(33) 



and similar expressions for Oq ,0?, , 9q ,63. In order to calculate the energy density and pressure of 
dark matter, let us consider a plane wave traveling in the x-direction 



(3 = B cos(ut — kx) ; (B, u>, k = const.) 
Making use of (pj) we get from (EG 



5 — k ~\~ k 
& 



In the quantum-mechanical representation one writes 

e = JvjJ , p = hk . 



(34) 
(35) 
(36) 



where e, and p are the energy and momentum of the particle, so that with ( |24| ) relation ( pq ) may be 
written as 



2 2 2 1 2 4 

s = c p +m c 



For the plane wave given by (|4j) we obtain from ( |3l| ) 



e° = (1/2) p 2 

x-wave 



U , 2 


(fr_ 




-5T + k 2 + 
c l 


yR 2 





(37) 



(38) 
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To describe an ensemble of particles having spatial isotropy, we can take a combination of six plane 
waves: the wave ([54]) moving in the positive direction allong the x-axis, one moving in the opposite 
direction, two waves moving along the y-axis in both directions, and two moving along the z-axis. If we 
take the time average, we obtain for the combined system of waves 



O = (3/2)B 2 _ + *» + _ + _ . (39) 




Making use of rt32) , we obtain for the same system of six plane waves 



©0 = 00 = 00 = . (40) 



Finally, making use of (18), (|22|), and pa), we get from equation (R3j 



e\ = el = el = -le° . (4i) 



The expression (39 ) is actually the energy density of our dark matter p D . One can consider it as consisting 
of two parts, the first being the contribution p B of scalar bosons, the second, p G , representing a global 
effect. 

P B = (3/2) B 2 + k 2 ^j , p G = (3/2) B 2 + ±^ , p D = e" = Pb + p G . (42) 

Up to here a monoenergetic ensemble of bosons is discussed. According to ( |25|) t he boson has zero 
mass in the radiation era, whereas in the dust period it is extremely light (cf. (|29|) , and J3C|)). These 
particles are discribed by the Bose-Einstein statistics, and there is no interaction with ordinary matter, 
so that the total number of bosons in the universe is conserved. In order to get a realistic scenario, one 
can turn to a gas of scalar bosons in thermal equilibrium, and fit it into the cosmological model. This 
procedure was already carried out in (Israelit & Rosen 1992, 1994) for a gas of weylons. 



5 A global approach. 

Above a F-R-W universe containing small regions of local inhomogeneities was considered. In the present 
section we consider a completely homogeneous and isotropic universe, with (3 depending only on the 
cosmic time t. With the metric given by flTI]), we can rewrite eq. (|l4j) as 

(/3i? 3 ) ^-^{p- iP )R^ . (43) 

Further, making use of (|l7|), and (|43|), we obtain from ( |l0| ) the following non- vanishing components of 
the energy-momentum density tensor of the global scalar field 

ei = ei = e| = ~e8 . (45) 

One can consider f3 as a function of the cosmic scale factor R, so that 

$ = R0,r with /3,r = ^ , (46) 

and ([U) may be rewritten as 

P D =e^ = ^(/3 + i?/3, fl ) 2 + |^ . (47) 
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It is remarkable that in this approach, like in the previous one (cf.(|3^)), the energy density p D of the 
scalar field is positive, so that one obtains a dark matter effect. 
Making use of (Q) we can rewrite eq. ( |43| ) as follows 

4"7T 

R 3 R 2 (3 :R , R + R 2 (RR + 3R 2 )P, R + —R 3 (p-3P)/3 = . (48) 

For a moment, let us go back to the relations (|l2|), and (|l5|). Making use of ( [Tt| ) we obtain for the energy 
density of dark matter 

fp D = ^6° = -g ; (23 = const.). (49) 
Let us consider the dust-dominated period. Here P <C p, and the energy relation ( pcif ) gives 



Substituting (|l8|), and (|l^ ) into (|4q), and taking into account (^), (p0|), we obtain 

[MR 2 + DR- R 3 ] [3 tRjR + 
and turning to the dimensionless variable 



^MR + 2D- 3R 2 



P,R = -\Mp , (51) 
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R 

R N 

with 2?at being the present value of the cosmic scale factor, we get 



[MR N z 2 + Dz-R 2 N z 3 ] (3,„ 



(52) 



-MR N z + 2D-3R z N z'\f3 tZ + -MR N P = . (53) 



The constants M, D, and R N apearing in ( |53| ) may be expressed in terms of two observable cosmological 
quantities, the present value of the Hubble "constant" H N , and that of the energy density of ordinary 
(luminous) matter p N . It is convenient to introduce a parameter r\ , given by 

= ■ (54) 

A second helpful parameter is the ratio 

x= P 1 ^ L ^ l+ Po_ (55) 
P P 

so that for the present state of the universe we have 

PtotaLN = PN + Pd.N = X N P N ■ (56) 

On the other hand, eq. (jlj) at present may be written as 

1 87T 

H » + W = Y {Pn +/Vn) ■ (57) 

With (13), and (||) this takes on the form 

Ptata, = P N + PD.N = PC ^ + ^\ ( 58 ) 

Finally, making in (||) use of @, ©, ©, and (||), we obtain 

(l + V 2 )R 2 N =MR N +D . (59) 

In a similar way we obtain from ( |56| ) 

X^M^ =MR N +D . (60) 
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From 



and 



and (pO) one obtains 



M 



1 + V 2 



D = 



1 



(i + v 2 )Rl 



Inserting ([5lj), and 



into eq 

3 



t) we get: 



z 



3z 2 



1±V1 z 2 

5(l+r, 2 ) 
2 X N 



(1 + rj 2 ) Z 



2(Xv-l) 



(1 + ry 2 ) 



(61) 
(62) 

(63) 



2x f 



/3 = 



For given values of parameters ry, and \ N one can obtain the dark matter function P(z) from equation 
(|63|). It is believed (Tyson 1992), (Blome, Priester,& Hoell 1995) that at present the total density of 
matter exceeds that of luminous by a factor of ~ 10. This justifies taking 



X N = io 



(64) 



From (|5J), and ( p8| ) we see that large values of rj describe a nearly flat universe. Below will be considered 
models with r\ = 1, 4, 20, combined with x N = 10. 

Let us take the Einstein equation (Eq). Making use of (|49j), and (pQ), we can rewrite it as 



i? 2 



i 

r 2 



M 



R* 



A closed universe achieves its maximum radius at R = 0, so that from 
Further, substituting ( |49| ) into(^), we obtain 

^ a r/9+Bfl ^ 4. ^ - 3 D 



(65) 

we can calculate R ma x- 
(66) 



Making use of the values of R m ax obtained from (pq), one can calculate /3 max = /3(R max ), the value at 
the state of maximum expansion. Finally, taking into account (^), and (|50|), as well as ([H]), and $3% ) 
we have for the ratio ( |55"| ) during the dust-dominated period 



so that 



x(R) = i 



X, 



D 
~MR 



1 



= 1 



(Xjv - 1)^ 

i? 



DR. 



Rn 



(67) 



(68) 



Taking the values of Rmax obtained from eq. (p5|), we can from fl68|) calculate the values of Xmax- For 
a moment let us go back to (^lj) . At the state of maximum expansion the first bracket in that equation 
vanishes (cf.j65|)). Thus, making use of (61), and (|62|), and of the above mentioned values of f3 max , one 
can calculate the value of the derivative (0,R) mnm , for the moment when R = R m ax- By the procedure 
described above, we get from equations (65), (^), (|68|), for chosen values of x N an d r], the results 
listed below: 



Xn = 10; 7] = 1; = 1.4i? N ; (5 max = 0.45; Xm„ = 7.4; (/?,fl) ma:c = 2.3 x lO" 2 ^ (69) 

Xjv = 10; r, = 4; = 4.8i? N ; f3 max = 0.39; Xm « = 2.9; G9,fl) mo!B = 4.3 x lO" 4 ^ 

X N = 10; V = 20; = 48i?„; /3 ma:c = 0.19; Xma* = 1.2; (/?,*)_„ - 1.2 X 10- 4 if N 



From (58) one can readily see that the first line in (|69j) belongs to a well closed universe, whereas the 
third line describes an almost flat model. It must be noted that we have made use of two equations, 
©, and ©, both for the same function f3(R). However differentiating ( |66] ) we can verify that these 
two equations are coexistent. Now, if one wants, he can obtain @(R) from the linear differential equation 
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( pl| ) (or alternatively from (|63|)) for a chosen pair of parameters x N , V- l n this procedure the above given 
values of (3 max ee f3(R max ), and (f3,R) max , serve as boundary conditions. 
Let us consider equation (|65|). Integrating it one obtains 

t + C = -(£> + MR- R 2 ) 1 ' 2 - — sin' 1 J^_jjL (70) 

with C being a constant. Inserting (|6l|), and (|62]), turning to z according to (|52|), and making use of 
one can rewrite (|7(i|) as follows 



<+c— 



2 lV2 u „2 



— {x N -l + z)-z 2 



sin -i i + ^-2 X „* . (n) 

2X„ (l + r^l + %^ 



With the help of (^Tj) one can estimate the age t n of the universe. For this purpose one neglects 
the duration of the radiation, and preradiation periods (cf. (Weinberg 1973), (Israelit, & Rosen 1989)). 
One also must take into account the fact that the radius of the universe at the beginning of the dust- 
dominated period Rq is much smaller then its present value R N , (Rq < 0.01R N ). Further, making use of 
the values of R ma x quoted in ([39]), one can also estimate the moment r maxi when the universe reaches 
its maximum radius. Straightforward calculations give the following results. 

X N =10; r, = 1; => t n = 0A2H- 1 ; W = 1.13H- 1 . (72) 
Xn =10; 7? = 4; => t n = 0.50F; 1 ; W = 14.21/- 1 . 
X N = 10: 77 = 20; => r N = 0.51H- 1 ; r max = 1447J- 1 . 



Let us take for the present value of the Hubble constant (cf. (f26|)) 

H N = 50 km s" 1 Mpc^ 1 = (2 x 10 10 y)" 1 = 0.5 x 10~ 28 cmT 1 . (73) 
Then we obtain from ( [72] ) for the models considered above 

0.84 x 10 28 cm < r N < 1.02 x 10 28 cm (74) 
and turning to time units we can rewrite this as 

8.9 x 10 9 y < r N < 10.2 x 10 9 y . (75) 



6 Discussion. 



In the present work dark matter is obtained from a conformally coupled scalar field. It is shown that this 
field may be treated as a source of scalar bosons and that it can also cause a global dark matter effect. 
In the process of creating dark matter the role of the conformally coupled scalar field is similar to that 
of the Dirac gauge function. However, whereas it was possible to choose the latter almost arbitrarily 
(cf. (Israelit, & Rosen 1995)), the behavior of the conformally coupled field dep ends on the density and 
pressure of ordinary (luminous) matter in the universe (cf. equations (|I^), (pl|), ([43])). 

Two aproaches are considered, a local, and a global one. In the first case we assumed that the globally 
homogeneous and isotropic universe contains small regions of local inhomogeneities, wherer the scalar 
function (3 may depend on spatial coordinates. This dependence leads to a Klein-Gordon equation for 
P (eq. (|23|)), so that from the quatum mechanical standpoint one obtains an ensemble of neutral and 
spinless bosons. During the radiation period these particles have zero mass, and in the dust-dominated 
universe they are extremely light, so that these particles constitute dark matter pervading all of the 
cosmic space. 

In the global approach (section 5.) the universe is completely homogeneous and isotropic, and the 
scalar function /3(t) generates a global dark matter effect. In order to specify cosmological models, two 
parameters are introduced, the first, ry characterizing the degree of flatness of the universe, and a second, 
X being the ratio of the total matter density to that of ordinary (luminous) matter. It is assumed that 
at present x N = 10, and for various values of 77 the present age of the universe, its maximum radius, as 
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well as T max , the moment when the expanding universe will achieve its maximum radius are calculated. 
From (§|), @-75) one can conclude that the model with x N — 10 and r\ — 20 (an almost flat universe) 
agrees with nowadays cosmological observations (Bolte, & Hogan 1995), (Chaboyer, Demarque, Kernan, 
& Krauss 1996). For this model we have: R m ax — 48i? N ; t n = 10.2 x 10 9 y; T max = 282r N . In the 
expanding universe, during the dust-dominated period, the ratio n (cf. (|67[)) decreases by a factor of 
about 1000, hence this dark matter form will not play an essential role in the late stages of the expansion 
phase, but it was important in the early stages of the dust universe. So it might have affected growing 
processes of galactic formations at the beginning of the dust era (Israelit, Rose, & Dehnen 1994, a, b). 

The conformally coupled dark matter considered in the present paper is to be regarded as a dark 
matter form connected closely to the Weyl geometry. Assuming that the space-time of the universe 
has Weylian properties (in addition to the well known Riemannian properties) we can derive several 
geometrically based dark matter forms (cf. (Israelit 1995)): a) the dark matter considered in this work, 
b) a gas of heavy weylons (Israelit, & Rosen 1992, 1994), c) a global dark matter effect from the Dirac 
gauge function (Israelit, & Rosen 1995). All these forms may be fitted into cosmological models, giving 
agreement with observation. 
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